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Abstract 

A comprehensive investigation is carried out on the class of energy based ensem- 
bles. The eight ensembles are divided into two main classes. In the isothermal class 
of ensembles the individual members are at the same temperature. A unified frame- 
work is evolved to describe the four isothermal ensembles. Such a description is 
provided both in the second and the third constraint formalisms. The isothermal- 
isobaric, grandcanonical and the generalized ensembles are illustrated through a 
study of the classical nonrelativistic and the extreme relativistic ideal gas models. 
In the adiabatic class of ensembles the individual members of the ensemble have the 
same value of the heat function and a unified formulation to described all the four 
ensembles is given. The nonrelativistic and the extreme relativistic ideal gases are 
studied in the isoenthalpic-isobaric ensemble, the adiabatic ensemble with number 
fluctuations, and, the adiabatic ensemble with number and particle fluctuations. 
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I Introduction 



Boltzmann-Gibbs extensive statistical mechanics was generalized by Tsallis through the 
use of a deformed logarithm function in the entropic expression [1]. The functional form 
of the entropy defined in terms of the deformed logarithm is 

S, = -kJ2pj\n,p,, gGM+, (1.1) 

i 

where k is the Boltzmann constant and pi is the probability associated with the i-th 
microstate. The positive real condition imposed on the nonextensivity parameter q ensures 
the stability of the Tsallis entropy [2\. The g-deformed logarithm made use of in fll.ip 
and its inverse the g-exponential are defined below 

lnga; = exp^(x) = [1 + (1 - g)a;]^, (1.2) 

where we have used the notation [x]+ = max {0,a;}. In the g — )■ 1 limit, (11. ip reduces 
to the extensive Boltzmann-Gibbs entropy. A detailed account on the properties of the 
g-entropy and its applications may be found in [3j. 

To describe a system in thermodynamic equilibrium with its surroundings, we need 
three macroscopic variables corresponding to thermal, mechanical and chemical equilib- 
rium. A system with fixed values of these parameters can exist in different microstates. 
A collection of systems existing in all the possible microstates corresponding to the same 
value of the macroscopic variables is called an ensemble. Each of these macroscopic param- 
eters can be chosen between an extensive and an intensive variable. Thus we have eight 
different ways of choosing the three variables and consequently eight different ensembles. 
The eight ensembles can be divided into two different classes namely the isothermal class 
and the adiabatic class. In the isothermal class of ensembles, the thermal equilibrium 
variable is the temperature and the individual members of the ensemble are at the same 
temperature. The canonical, isothermal-isobaric (HE], grandcanonical and the generalized 
ensemble [1|5] belong to the isothermal class. Each one of these isothermal ensembles is 
characterized by an heat function and a free energy. A general expression for the heat 
function and the generalized free energy are given below. 

= U + J2^iXe (1.3) 
JF = fZ-TS + ^x^X^ (1.4) 

where X is an extensive thermodynamic quantity whose intensive conjugate is x. The 
index i can take three different values of either 0, 1 or 2. In the i = case the internal 
energy and the Helmholtz free energy are the heat function and the free energy respec- 
tively. The ensemble pertaining to this situation is the canonical ensemble {N,V,T). 
Both the isothermal-isobaric {N,P,T) and the grandcanonical {fi,V,T) ensemble relate 
to the i = 1 case in (11.31) and (11.41) . The enthalpy H = U + PV and the Gibbs free 
energy G = U — TS + PV are the heat function and the free energy corresponding 



1 



to the isothermal-isobaric ensemble. For the grandcanonical ensemble the Hill energy 
L = U — fiN and the $ = f/ — TS — fiN are the heat function and the free energy respec- 
tively. The i = 2 case describes the completely open system {fi,P,T) and the ensemble 
is usually referred to as the generalized ensemble. The R function R = f/ + PV — fiN 
and S = U + PV — TS — fiN are the heat function and the free energy of the gener- 
alized ensemble. The adiabatic class of ensembles arise when the thermal equilibrium 
happens with respect to the heat function which is an extensive quantity. Based on the 
four different heat functions namely U,H,L and R we have four different ensembles. The 
well known microcanonical {N,V,U) ensemble. The isoenthalpic-isobaric {N,P,H) en- 
semble introduced in [7] and later extensively studied in references p|9] . More recently 
the adiabatic ensemble with number fluctuations {fi,V,L) and the adiabatic ensemble with 
number and volume fluctuations (/x,P,R) were introduced through references [lU] and [T^ 
respectively. All the four adiabatic ensembles were treated in a unified way in [12j. Of 
the eight ensembles, the microcanonical, the canonical and the grandcanonical ensem- 
bles have been extensively studied and are well known in statistical mechanics. But the 
other five ensembles are also of fundamental importance in studying physical problems 
whose experimental conditions emulate the conditions of these ensembles. For example 
the isothermal-isobaric and the isoenthalpic-isobaric ensembles are used to study fluids 
which are confined at constant pressure. A molecular dynamic simulation for these ensem- 
bles was evolved in [13]. The generalized ensemble has been studied in [H]. Applications 
of the generalized ensemble include lattice gas models, and a system of polymer molecule 
of monomers held together by weak forces such that the number A^ fluctuates [H]. 

An investigation on the ensemble formulation of nonextensive statistical mechanics 
based on g-entropy (11. ip was initiated in [1] through the study of two level system in the 
microcanonical and the canonical ensembles. A further detailed study of the canonical 
ensemble was carried out in references P^fTB] wherein the second and the third constraint 
models were introduced. In the second constraint method P3] we use the unnormalized 
g-expectation value which, for a given observable O reads: 



where pj is the ensemble probability of the microstate j in the second constraint. Though 
the thermodynamic Legendre structure is preserved in the second constraint, the formal- 
ism suffered from certain disadvantages. The unit operator does not preserve its norm 
for an arbitrary value of q and, the form of the energy conservation principle in the 
macroscopic and microscopic limits varied. In order to overcome the disadvantages in the 
second constraint, an appropriately normalized form of the g-expectation value known as 
the third constraint method was introduced in [TB]. The definition of the g-expectation 
in the third constraint formalism is 



Though the third constraint method cleared all the problems previously encountered 
in the second constraint method, the probabilities became implicit quantities rendering 




(1.5) 




(1.6) 
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the calculations difficult. But fortunately, the second and the third constraint models 
can be interrelated via a temperature transformation relation which helps us to override 
the difficulties caused by the implicitness in the probabilities and the thermodynamic 
variables. The various applications of the canonical formulation may be found in pTIITS] . A 
similar formulation of the grandcanonical ensemble was ffist done in the second constraint 
[T^ and later in the third constraint [20j. A detailed investigation of the microcanonical 
ensemble based on the g-entropy f ll.ip was carried out in [2T] and the classical ideal gas 
was discussed as an application. Thus in the nonextensive statistical mechanics, based on 
g-entropy only three ensembles, the canonical [T|15fl6] the microcanonical [TpT] and the 
grandcanonical p^|20] have been investigated so far. 




In the current work we present two unified frameworks one for the isothermal class 
of ensembles and the other for the adiabatic class of ensembles. The unified framework 
pertaining to the isothermal class of ensembles has been presented for both the second 
and the third constraint formalisms though the latter is the currently accepted formula- 
tion. Adequate explanation for providing a unified formulation for the second constraint 
formalism is discussed at the relevant places in the article. A generalized formulation of 
the temperature transformation relation has been constructed to interrelate the second 
and the third constraint formulations. An exact solution of the specific heats at con- 
stant pressure and at constant volume corresponding to both the classical nonrelativistic 
ideal gas and the extreme relativistic ideal gas has been obtained in the isothermal iso- 
baric ensemble. Contrarily obtaining an exact solution to these ideal gas models in the 
grandcanonical and the generalized ensemble appears to be difficult. We make use of the 
technique of disentangling of g-exponential into a series of ordinary exponentials [22] to 
construct a perturbative approach to the problem. Such an approach has already been 
used in the context of the canonical ensembles [T8|23] . In the case of grandcanonical en- 
semble, the perturbative terms up to second order in the expansion parameter (1 — q) has 
been retained. The perturbative terms have been displayed up to ffist order in the expan- 
sion parameter (1 — g) in the case of the generalized ensemble. The limitations imposed 
on the order of (1 — g) is only for the sake of simplicity and in practice the method could 
be extended to any arbitrary order in the perturbative parameter. A unified framework 
to describe all the four adiabatic ensembles in the nonextensive g-statistics has been put 
forward. A g-generalization of the equipartition theorem and the virial theorem have also 
been discussed. The specific heat at constant pressure and at constant volume correspond- 
ing to the nonrelativistic and extreme relativistic models of ideal gas have been exactly 
computed. Since a summation over the number of particles could not be evaluated in the 
(/i,y,L) and the (/i,P,R) adiabatic ensembles, the phase space volume and the equation of 
state are expressed as formal sums. Throughout the article we identify the temperature 
with the Lagrange multiplier corresponding to the constraint on the internal energy. This 
is because if we assume the notion of quasireversibility and simultaneously implement the 
ffist law of thermodynamics the Tsallis entropy becomes identical to the thermodynamic 
Clausius entropy ^24j only when the Lagrange multiplier /3 is associated with the inverse 
temperature. 

The plan of this article is as follows: The unified framework for the isothermal ensem- 
bles is presented in Section [Tll wherein the second and third constraint formulations are 
described as separate subsections. Using the classical nonrelativistic and extreme rela- 
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tivistic ideal gases the isothermal-isobaric ensemble, the grandcanonical ensemble and the 
generalized ensemble are illustrated in Sections IIIH IIVI and |Vl A generalized formulation 
of the adiabatic ensembles is given in Section I VII The ideal gas model in the isoenthalpic- 
isobaric ensemble have been worked out in Section IVIII In Sections IVIIII and IIXI the 
{fi,V,L) and the (/i,P,R) ensembles are discussed. We present our concluding remarks in 
Section |Xl 

II Isothermal ensemble 

The isothermal class of ensembles is used to study systems in which the thermal equili- 
bration is with respect to the temperature. A unified framework to study the isothermal 
class of ensembles in introduced in this section. This section is divided into two sub- 
sections. In the first subsection a unified formulation of the unnormalized g-expectation 
values or the second constraint formalism is presented. A similar unified method of the 
third constraint formalism is presented in the second subsection. Finally we demonstrate 
the interrelation between the two formalisms and derive a generalized expression for the 
temperature transformation relation. 

A. Second constraint 

The expectation values of energy, a general thermodynamic observable and the heat func- 
tion based on the unnormalized g-expectation value (11.51) is 

[/f (Xi,X2,/3) = 5^ (pS^^^(Xi,X2,/3))' e„ (2.1) 
xg(Xi,X2,/3) = 5^ (pS^^^(Xi,X2,/3))' X,. (2.2) 
^f)(Xi,X2,/3) = 5^ (pS^^^(Xi,X2,/3))' (e, + 5^x,X,) (2.3) 

(2) 

where Pj ^^^j is the probability of finding a particle in a particular microstate and X is 
an extensive thermodynamic quantity whose intensive counterpart is x. The chemical 
and the mechanical equilibrium variables are denoted by Xi and X2. A brief introduction 
of the variable i and the method to deduce the various ensembles corresponding to its 
different values is given in Section [H 

Optimization of the entropy (11. ip using the normalization condition, the internal en- 
ergy (12. ip and the other extensive thermodynamic variables (12. 2p leads to the distribution 
given below 

pS,.> (^1' ^2, /3) = \ exp, ( - /3 (^^ + E ^^^^)) ' (2-4) 



4 



(2) 

where Zq the partition function in the second constraint is 



(2.5) 



As previously stated in section HI we associate the temperature with the Lagrange multi- 
pher corresponding to the internal energy i.e., /3 = 1/kT. 

The second constraint heat function of a given ensemble {Xi, X2, (3) can be computed 
from the partition function using 



iof (Xi,X2,/3) = -^ln,Zf (Xi,X2,/3). 

In terms of the partition function the extensive variable in the second constraint is 

1-9 



/3 dxe 



ln,Zf (Xi,X2,/3). 



(2.6) 



(2.7) 



The knowledge of the heat function and the extensive thermodynamic variables of a given 
ensemble enables us to find the internal energy using Legendre transformation. 

The thermodynamic potential corresponding to an ensemble in the second constraint 
formulation is the generalized free energy 

J-f = -fcTln,Zf (Xi,X2,/5), (2.8) 

from which the entropy and the other extensive variables are computed as follows: 



S = 



dT 



^(2) _ Oj-g 



dxe 



(2.9) 



The thermodynamic expression for the specific heats corresponding to the different con- 
ditions is 



^(2)1 ^ 
^1 \v QT 



•{1} 



V 



rji y 



dT 



(2.10) 



B. Third constraint formalism 

To overcome the failures of the second constraint we use the normalized g-expectation 
values in the third constraint formalism. The definition of the internal energy, an arbitrary 
extensive variable and the heat function in this formalism are 



[/f (Xi,X2,/3) 



^tgi'^iy -^2, /3) 



f|(=^)(Xi,X2,/3) 



1 


c(3)(Xi,X2 


P) 


1 




c(3)(Xi,X2 


P) 


1 




c(3)(Xi,X2 





j2 {p\%,,,i^i,^2,p)y Q, (2.11) 

(pS,,>(^i,:£2,/3))' Xe, (2.12) 
J2 (pS,,j (-^1, ^2, /?)) + ^iX,) ,(2.13) 
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(3) 

where Pl-x^iy is the probabihty of finding the particle in a particle microstate. The factor 
c*-^-* is the sum of g-weights in the third constraint and is defined as 

c^'\X„X„P,)= J2 (pS^.,(^1'^2,/3))'. (2.14) 



Along with the normalization condition, the definitions corresponding to the internal 
energy fl2.1ip and the arbitrary extensive variable (12.121) are used in the optimization of 
the entropy (11.11) . The probability distribution obtained through this procedure is 

^' ^ if(^k« ( - ^^d^) + ^ - ^'') ) • 

(2.15) 

where S^q is the heat function in the third constraint and f3 = 1/kT. The generalized 
partition function Zq^^ in (I2.15P is as defined below: 



16) 



The generalized partition function and the sum of the g-weights are related via the ex- 
pression 

(zf)(/3,Xi,X2))'"'' = c(3)(/3,Xi,X2), (2.17) 

which holds for all the isothermal ensembles and plays a crucial role in solving the implicit 
equations in the third constraint. 

The heat function in the third constraint satisfies the differential equation 

/3 (/3, ^1, ^2) = ^ In, (/3, ^2). (2.18) 

From the knowledge of the heat function i^g , and the generalized partition function Zq 
the expectation value of the extensive thermodynamic variable is found using 

xg(/3,Xi,X2) = ^ (^^f (/5,Xi,X2) - ^ln,Jf (/3,Xi,X2)^ . (2.19) 

The internal energy can be obtained from the heat function and the extensive thermody- 
namic quantities using Legendre transformation. Making use of the property (12.171) . the 
nonextensive entropy (II. ip expressed in terms of the generalized partition function (I2.16P 
is 

S', = fcln,Zf (/3,Xi,X2). (2.20) 

Substituting the entropy (I2.20p in (11.41) we arrive at the free energy in the third constraint 
as 

= - iln, Jf)(/3,Xi,X2) + Y.x,xfl. (2.21) 
P {1} 
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In terms of the generalized free energy defined in (12.21 1) the heat function and the q- 
expectation value of the extensive thermodynamic variable reads: 

The heat capacities at constant pressure and at constant volume in the third constraint 
are 



(3) 



V dT ' 

d^f dSq ^d^J'c'^^ 



(2.23) 



T — - = -T —. (2.24) 



The third constraint approach is the currently accepted formulation of nonextensive 
statistical mechanics and problems which do not involve simulations can be solved directly 
in this approach. In the case of Montecarlo and moleculardynamics simulations direct use 
of the third constraint is difficult. To overcome this, the calculations are initially per- 
formed in the second constraint and the results are then transferred to the corresponding 
quantities in the third constraint. The transformation relation for the canonical ensem- 
ble was introduced in [16]. Adopting a similar procedure for the unified framework of 
isothermal ensembles, the probabilities in the second and the third constraint formalisms 
are related as follows 



P^X^iP, Xi, X,) = pS^,j(/3', Xi, X,) = exp^ ( - /3'(e,, + $^x,X,)) (zf (/?', Xi, X^))"' , 

(2.25) 

where, the pseudoinverse temperature f3' is related to (3 as 

= ^ . (2.26) 

c(3) (/3, Xi, X2) + (1 - g) /3 Sj^i^ [P, Xi, X2) 

The relation in (I2.26P is inverted to obtain /?' in terms of /3 

B = B' (2 27) 

1 _ n - «^ (/3'.^i.^2) 
J- i-L c(2){/3',xi,x2) 

After obtaining the transformation pertinent to a given system, there are two independent 
methods through which we can get the other thermodynamic variables. 

(i) Probability transformation technique: In this method, introduced in Ref. [H], the 
sum of g-weights is transformed from the second constraint to the third constraint using 
the temperature transformation relation (12.271) . Using (I2.17P the generalized partition 
function is obtained from the sum of g-weights from which we can calculate the specific 
heat directly. 

(ii) Observable transformation technique: Any dynamical observable in the second 
constraint can be transformed to the corresponding one in the third constraint [23] through 
the relation given below 

(2.28) 
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All the relevant thermodynamic quantities of a given ensemble in the third constraint can 
be computed from the second constraint using the two methods described above. 

The thermodynamic relations for the canonical ensemble can be obtained from the 
unified framework by observing that the internal energy and the Helmholtz free energy 
are the heat function and the free energy respectively. An elaborate exposition on the 
canonical ensemble may be found in P^fTB] and since this ensemble has been investigated 
in detail we do not intend to discuss it anymore. 

Ill Isothermal-Isobaric ensemble 

The isothermal-isobaric ensemble is used to study systems which attain equilibrium by 
exchanging both internal energy and volume with the heat bath. The enthalpy and the 
Gibbs free energy plays the role of heat function and free energy respectively. Using the 
unified framework described in Section HIl the nonrelativisitic and the extreme relativistic 
ideal gases are investigated. 

A. Nonrelativistic Ideal gas 

The Hamiltonian of a classical ideal gas in the non-relativistic regime with particles pos- 
sessing momenta Pj (i = 1, A^) in D dimensions is 

2 

^nr{p)=J2^^ P.= |P.|- (3.1) 

i 

The generalized partition function corresponding to the third constraint in the (A^, P, T) 
ensemble is 

N 

X "-p, (-|y (e ll - uf' + Py - Py^') ) ■ (3-2) 

Carrying out the phase space integration in (13. 2p . we arrive at 

^'"<^- ^' ^ ^ 70k) ' " ' 

P,3) 

where Ai = (^p^)^''^ and A = 1 + (1 — q)DN/2. In order to express the thermodynamic 
variables in a simplified manner, we define the following quantities 

A, = 1 + (1 - = (1 - (1 - fi. = A,£«, (3.4) 

where k = 1 (2) for a non-relativistic (extreme relativistic) case and 'Di = + A^ 
{T>2 = DN + A^). Integrating over volume in (13.31) we arrive at the generalized partition 
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function 



/3 



Ai-(l-g) 



(3.5) 



It is interesting to note that (13.51) has the dimensions of volume. This is due to the 
overcounting of volume eigenstates in systems in which the volume is a continuous variable. 
To overcome this we use the shell particle method of counting of volume states [25|26] . 
In the shell particle method of counting we take into account only the minimum volume 
needed to confine a particular configuration. When there are many equivalent ways of 
choosing a minimum volume we treat all of them as the same volume eigenstate. The 
minimum volume is fixed by enforcing a condition where we require that at least one 
particle lies on the boundary of the system. Carrying out the volume integral of the 
generalized partition function using the shell particle counting technique to reject the 
redundant volume states we arrive at 

(/3,P) = M- (exp, (|^(f/f + Prf )))'^ (^^y^ (i)". (3.6) 

The above expression of the generalized partition function (13. 6p in contrast to (13.51) . is 
dimensionless making it consistent with the basic principles of statistical mechanics. The 
ratio of gamma functions Qi made use of in (13. 6p is 



(i-gF^r(fEf + Di)- 

The heat function of the isothermal-isobaric ensemble given by the basic relation (I2.13P 
is computed and the result is 

-f^^.^e. u(4.H'AY'(^r('^r'. (3.8) 



The ratio between the (13. 6p and (13.81) is computed 



1/(3) = 2)^ (3.9) 



(3 



Using (13. 9 p in conjunction with (I2.17p . in (13. 6 p the explicit form of the generalized partition 
function is 



Yi^^ = {M- g^f {e.p^mr [jj (Ij . (3.10) 

Similarly the explicit expression of the enthalpy computed through the use of (13. Sp . (I3.10p 
and fl2T71) is 
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where the constant factor H„r used in (I3.1ip is 



(3.12) 



The entropy of the system in the third constraint written down from the generahzed 
partition function using the relation fl2.20p is 



k 



1-q 



(l-<?)Di£i 



(l-g)Aff 



- 1 



(3.13) 



Making use of the thermodynamic Legendre structure fl2.2ip the Gibbs free energy of the 
classical ideal gas is obtained: 



(l-g)/3 



((l-g)2)i-l) 



(l-<?)Ar£i 



- 1 



(3.14) 



The average volume of the classical nonrelativistic ideal gas in the third constraint can 
be found via the expression fl2.19p using the explicit expression of the generalized partition 
function fl3.10p and the enthalpy f l3.1ip . The calculated expression of the average volume 
reads: 



y(3) 



(l-<?)Di£i 



(l-g)Af£i 



(3.15) 



The average volume which was computed using enthalpy can also be calculated through 
two other different methods. The first method is based on a specific form of f l2.22p 
applicable to the isothermal-isobaric ensemble in which the average volume is related 
to the Gibbs free energy. The second method makes use of the defining relation of the 
average volume fl2.12p to obtain an implicit expression. Using the implicit expression in 
conjunction with (13. lip . (I3.10p and (I2.17P we can obtain the final explicit form. The results 
obtained through both these methods are in agreement with that obtained in (13.151) . 
Rewriting the expression for the average volume (I3.15P we arrive at the nonextensive 
form of the equation of state corresponding to the isothermal-isobaric ensemble 



py(3) 
9 



(l-g)Di£i 



{l~g)Af£i 



(3.16) 



From the knowledge of the enthalpy and the average volume, the internal energy evaluated 
using the Legendre structure is 



DN _ 



(l-q)Nl 



(3.17) 



The defining relation of the internal energy (12. lip has been used to verify (I3.17P in a 
manner analogous to the calculation of the average volume. 

The specific heat at constant pressure can be found from either the enthalpy (13. lip or 
entropy (I3.13P or the Gibbs free energy (I3.14p using (I2.24p 



DN 



(l-<j)X)i£i 



{l-q)Ni 



(3.18) 



10 



From the internal energy fl3.17p the specific heat at constant volume is evaluated through 
the use of fl2:23|) 



Cf|,, = ^*=„, (i) (i) . (3.19) 

The difference between the specific heat at constant pressure flS.lSp and the specific heat 
at constant volume (13.191) is 

CiX-C!,\ = Nk [i^^ (^ij , (3.20) 

which is the generalized form of Meyer's relation and in the extensive limit the Ihs goes 
as Nk. It is interesting to note that the factor 7 expressing the ratio between the two 
specific heats 

^1 \v 

is independent of the nonextensive parameter q and the number of particles N, though 
the specific heats (13.181) and (I3.19P depend on these parameters explicitly. 

The thermodynamic quantities calculated in this section can also be evaluated in a dif- 
ferent method where we compute the thermodynamic quantities in the second constraint 
framework. The ensemble probabilities in the second and the third constraints are related 
through a fictitious temperature /3' as given in equation (12.250 . This fictitious temperature 
P' in the second constraint is related to the original temperature /3 in the third constraint 
through a transformation. The temperature transformation relation corresponding to the 
classical ideal gas in the isothermal-isobaric ensemble is 

j; = {M- g,f-^^'^ Al'^ [^j . (3.22) 

Using the above transformation the third constraint thermodynamic quantities can be 
computed using the probability transformation technique or the observable transformation 
technique. The thermodynamic quantities and the equation of state obtained through this 
transformation method are identical to the corresponding quantities computed directly in 
the third constraint. 

B. Extreme relativistic ideal gas 

The Hamiltonian of the extreme relativistic classical ideal gas in D dimensions is 

7ier{p) = C ^Pi, Pi = \pi\, (3.23) 

i 

where the factor c is the velocity of light. The nonextensive study of the Hamiltonian 
(I3.23P is carried out in the isothermal-isobaric ensemble. Since the calculational proce- 
dure is identical to the nonrelativistic case we present only the core steps leading to the 
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calculation of the specific heats. The generalized partition function calculated using the 
basic definition is 

= (^) ) . (3.24) 

The notations made use of in (13.241) are 

^ - (ch)or(D/2) ' - (1 _ r (f5f + s,-^ ' ' 



Equation (I3.24p is an implicit equation since it is a function of the sum of g-weights 
which is related to the generalized partition function through the relation (12.171) . The 
nonextensive enthalpy of the extreme relativistic ideal gas is 

,,3,.,,^,Je.p,(A,,3AyVn7a-- (3.26) 



where the factors A2, £2 and fi2 used in this section are described in (13. 4p and the discus- 
sions following it. The ratio between (I3.24p and (I3.26P is 

r(3) 

i/f = 2)2 — . (3.27) 

Using fl3.27p in fl3.24p we obtain the explicit expression corresponding to the generalized 
partition function 

= (A^g,)'^ (exp,(D.))- (^ij . (3.28) 

Similarly the calculated value of the explicit form of the enthalpy is 

^ / 1 \ a~q)NZ2 / -I \ (l-g)Z)2£2 

= ^ -e. ( 4 ) I ^ ) , (3.29) 



where the factor Sgr is defined as 

Sg, = (A^^?2)^'"'^''' . (3.30) 

The average volume of the extreme relativistic ideal gas is evaluated through the use of 
(16. 7p . The equation of state obtained by rewriting the expression for the average volume 
is 

(^) . (3.31) 

To find the internal energy we use the Legendre transform (II. 3p and the expressions for 
the enthalpy and the average volume. The computed expression of the internal energy 
rGcicis ' 

= — =.4^) [p) ■ (3.32) 
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The specific heat capacity of the extreme relativistic ideal gas at constant volume and 
at constant pressure can be established from the final explicit expressions of the internal 
energy and the enthalpy respectively. We exhibit the results below 

^ (l-q)2)2£2 / ^ \ (l-q)NS2 



CiX = DNk S,, 2, (^-j (^-J , (3.33) 

/ . N (l-q)S)2£2 / . \ {l'q)N£2 

CfL = {DN + N)k E^rS., {-) - . (3.34) 



■q IP - y^j yp 

The ratio between the heat capacity at constant pressure and at constant volume 

c^fL . 1 



1 + -, (3.35) 



^1 \v 

is found to be independent of both the nonextensive parameter q and the number of 
particles A^. All the results in this section have been calculated through two methods. 
In the first method, the third constraint generalized partition function, sum of g-weights 
and the enthalpy are computed. Since these quantities are implicit in nature, they are 
systematically solved to obtain the explicit forms. The second method is based on the 
interrelation between the second and the third constraint formalisms. The thermody- 
namic quantities are first obtained in the second constraint formalism and they are later 
transformed to the corresponding quantities in the third constraint through the use of 
the temperature transformation relation. Both the probability transformation and the 
observable transformation techniques are used in this objective. The results obtained 
through these methods are in precise agreement with each other. The proper counting 
of volume states advocated in references (25126] has been used in both the methods. In 
the extensive g — )• 1 limit, the standard physical quantities corresponding to the classical 
Boltzmann-Gibbs statistics are recovered. 

Another interesting model is the Tonks gas [27] which is a one dimensional model of 
hard rod gas. The Hamiltonian describing the dynamics of a Tonks gas is 

2 



^*5(P) = E^ + ^W' P^ = \vl (3.36) 

i 

where the potential V{r) is 

V{r) = x>a (3.37) 
= oo X < a. 



An exact solution of the Tonks gas in the free length ensemble can be read from the 
results of the nonrelativistic ideal gas. In this discussion the term free length refers to 
the actual volume available for a particle Lf -q = L — {N — l)a. The equation of state of 
the Tonks gas in the free length ensemble can be read from the expression corresponding 
to the classical nonrelativistic ideal gas by observing that the number of dimensions is 1 
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and replacing the volume by the free length. The equation of state of the Tonks gas in 
the free length ensemble is 

AT . s /I \ ^-21:1 / 1 \ ^-^1:1 

= f (Ai:,)- (i) {^) , (3.38) 

where the functions Qi.i and Ai^ are special cases of (13. 7p . obtained by substituting D = 1 
and All is the 1 dimensional analog of A^. A direct computation of the Tonks gas is used 
to verify the expression for equation of state. The correct counting of volume states has 
been applied and the Boltzmann-Gibbs results are recovered in the g — 1 limit. 



IV Grand Canonical ensemble 

A system which exchanges both the internal energy and the number of particles with its 
surroundings is described using grand canonical ensemble. The Hill energy L = U — fiN 
is the heat function and the free energy is $ = ?7 — TS — fiN. In this section we extend 
the perturbative formulation developed in [18] and [23] to study the nonrelativistic and 
extreme relativistic ideal gas in grandcanonical ensemble. 



A. Nonrelativistic classical ideal gas 

A perturbative study of the nonrelativistic classical ideal gas described by the Hamil- 
toninan (13. ip is carried out in the grand canonical ensemble. The generalized partition 
function in the third constraint is 

Zf>(,, VJ) = f; ^ / d°"r d-p exp, U (e - - Lf)] . (4.1) 



N=0 



where ^ = and c'-^-' is the sum of q weights corresponding to the grand canonical 
ensemble in the third constraint. The Hill function Lg = Uq — ^Nq is the heat 
function corresponding to the grand canonical ensemble where, Uq and Nq are the 
internal energy and the average number of particles respectively in the third constraint 
and /i is the chemical potential of the system. 

The observation in [22], that the q exponential can be written as an infinite multi- 
plicative series of ordinary exponential 



exp„ {X) 



n-p((-i^^^^^4^)' (4-2) 



k=l 



which in turn can be expanded up to any order in (1 — g), is used to construct a perturba- 
tive procedure for the generalized partition function. Assuming the generalized partition 
function (14. ip to be well defined in the region g = 1, we replace the g-exponential by the 
ordinary exponential. The resultant expression for the generalized partition function is: 

)(/., V, /3) = V,{d,) exp(£) / d^^r d^> exp - ^( £- - ^At) . 

(4.3) 
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where C = l3Lq and the operator valued series 'Di{dp) in (14 .Sp is 

m.) = 1 - ^ + ^ (4'' + 1 i';) + ■ ■ • . = -3"^. (4.4) 

In the expression fl4.3p . we assume that the summation over the number of particles is a 
convergent, atleast in the region where the nonextensive parameter is close to 1. Carrying 
out the phase space integration and the summation over the number of particles in fl4.3p . 
the generalized partition function in the third constraint can be expressed as 

Zf\^^,V,l3)=V^{dp) exp(£) ZbgCP). (4-5) 

where Zbg{(^) = exp (exp(/3/i) ^ is the Boltzmann-Gibbs partition function. 

The perturbative series of the generalized partition function (14. 5p calculated upto second 
order in the expansion parameter (1 — g) is 

) = exp(£) ZbgCP) (1 + (1 - g) ^1 + (1 - g)' ^2 + ■ ■ ■ ), (4.6) 
where the coeffecients of the perturbative series (14. 61) are listed below 

+4(8 {hf - ^D') «(^) + ^ (4 - D')) + ^ {NbgCP)? (48 C {a{j3)f 

+8 £ (3 £ + 2)(a(/3))2 - 12D a(/3) + 8 (4 - 7D)0ij)^ + 42 0ij,)^ 

-2D^ (9 + 7D)^ + D\3 + 5D + 7D^)) + ^ NbgW) f 12 (a(^))' 

+48 £ {a0))^ + (72 £ - 48) £ {a{/3)y + (48 + 72 - 12 D^)a0) 

+64 0i^f - 12D Cp^f + 2AD /3fiC- D{A + 9D- D^)) , (4.7) 



The factor a0) = f3 — ^ and the average number in the standard Boltzmann-Gibbs 
statistics is Nbg0) = exp(/3yu) V Ad (5'^^"^ are used in (14. 7p . We notice that the general- 
ized partition function described above is an implicit quantity and depends on the heat 
function \Sq^ and the sum of g-weights c*^'^^ which is in turn related to the generalized 
partition function via the relation (I2.17p . 

The integral form of the sum of g-weights written down from the defining relation 
flXnp is 



^^'H/^, V^, /3) = T-T^-^ X (4.^ 



X 



tj^j d""" d™P (exp, ( - ^(5: ii - - Lf.)))' 
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A perturbative computation of the sum of g-weights is evolved in the following manner, 
first the g-exponential is expressed in terms of the ordinary exponential 

= -^—^ ^ ©2(^^/3) exp(£) X (4.9) 

(/i,V^,/3) 



where the operator series V2{dis) appearing in fl4.9p is 

V,(d,) = 1 - (1 - ,) (4" + i df) + (1 - ,f (df + I df + \ 4"') + . . . . (4.10) 

In the second step the phase space integral and the summation over the number of particles 
is carried out in fl4.9p and the resultant expression is 

^S^\^i,V,P) = ^—^ ^ V.idp) exp(£) Zbg0). (4.11) 



In the final step, the rhs in f l4.11l) is computed and the series corresponding to the general- 
ized partition function (14. 6p is substituted to obtain the perturbative series corresponding 
to the sum of g-weights 

c = l + (l-g)q3i + (l-g)2^2 + ---, (4.12) 
where the perturbative coefficients and ^2 are 

= -NBG0){a0)-l) 

^2 = -lcic-i) + ^iNBGmH^H~p))' + ic^m' + {^-D)c^m (4.13) 

+ \ NbgCP) (2 W)f + 2 (2 £ + 1) («(^))2 + 6D «(^) - (2 £ - 1) (4 - /}) + 4 C) 

The integral form of the grandcanonical heat function, the Hill energy Lq in the third 
constraint picture is 

'''''' - (ifi;^^^/^''"^^''"^^?--^'^) ' 

X (-P,(-KE|i-^'V-L<3.)) 

Using a perturbative procedure akin to the one employed in the sum of g-weights the Hill 
energy in terms of the derivative series f l4.10p is 

Lf(/x,r,/3) = j-—^ ^ V,{dp) {^-1^ exp(£) NbgCP) ZbgCP). (4.15) 
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Calculating the derivatives in the above expression and substituting the Boltzmann-Gibbs 
partition function, the perturbative series of the Hill energy upto (1 — g)^ order is 

Lf = iBaCm + (1 - g) Li + (1 - g)H2 + ■ ■ ■ ). (4.16) 
where Lbg(/3) = (yYp ~^ ^^ Nbg{P) and the coefficients of perturbation are listed below 

o 

+2D - 8 - ^ NBcii^) {2-D + 4D^f3fi-8 {13 ^if - 6 {a{j3)f + (3D 

-2D^ + 8 (1 - /}) - 8 (^/i)2) C)-l {NbgCP)? i^D^ - 16D P + 20/3^ 

8 

-16(«(^))3). (4.17) 

From the perturbative solutions we notice that the generalized ensemble (14. 6p . the sum of 
q- weights (14.121) and the Hill energy (I4.16P are coupled equations. The explicit expression 
corresponding to these quantities are obtained by employing a recursive procedure to 
solve these equations. The important characteristic feature that all these equations are 
uncoupled at g = 1 enables us to obtain a solution. Through the use of the recursive 
procedure, an explicit form of the generalized partition function is obtained up to second 
order in (1 — g) and reads: 

(/3) = ZBcim + (1 - ZiW) + (1 - + ■■■)• (4.18) 
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where the coeffecients of perturbation are 

MP) = -^{4D + AD^ + D^-(AD-2D^)I3ii+(8-AD)(Pii)^ + 8(/3ii)^)NbgW) 

+ ^{4D + AD^ + {8D + AD^)l3f, + 8(1 + D) {Pfxf - 8(/3/x)3) (iVBG(/3))', 

Z2{P) = (32L> + 88L>^ + 28D^ + 2L>^ + 3D^ - {32D + 168L'^ + 40D^ 

768 

-18D^) /3fi + {96D + 96^2 - 24D^) {13 - (512 - 416D + 240D-'^) {(5^if 

-(544 - 432L») (/3/x)^ - 96 (/3//)^) Nbg{P) + -^{mO^ - 32L>^ + 8^^ + 

512 

-(160D^ - 432D=^ - 8D^ - AD^) ^/i + {576D - 656D'^ - IdOD^ 
+12D^) (pii)^ - (1024 + 640D - 704D^ + 32D^) {pij,f - (704 + 1024D 

-A8D') (/3/x)^ + (384 - 64D) (/3/i)^ + 64 (/^/x)^) (iVsG(/9))' - ^ (32L> 

+48^2 + 16D^ + 2D^ + - {32D - 32D^ - 9QD^ - 2AD^ + QD^) (5^ 
+ (192 - 64D - 2081)2 - 1521)3 - 20/)^) (/^/i)^ - (320 - 128D - 368D^ 
+48D^) {Puf - (128 + 448D - 8QD'^){P^Y + (192 - 96D) (/3/i)^ 

-64 {I3nf) {NBG{P)f + 7^ {IQD'^ + 32D3 + 24D^ + 8D^ + - {64D'^ 

128 

+96D=^ + 48L>^ + 8D^) + (64D + 1920"^ + 144^^ + 32D^) {^/if - {192D 
+256D^ + 80D^) {P^if + (64 + 256D + 128D'^) {^^if - 128(1 + D) {^nf 
+U{^lif){NBGm^- (4.19) 

A similar application of the recursive procedure enables us to calculate the explicit expres- 
sion of the Hill energy upto second order in (1 — g). The perturbative series corresponding 
to the Hill energy is 

)(/3) = Ug(/5) (1 + (1 - q) Li(/3) + (1 - qf U{P) + •••), (4-20) 
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and the coefficients of perturbation are listed below 

Li(/3) = -( ^ ( + 2D^ + 12D f3fi + AD f3fia{f3) + 16 {(3fiy + 8{(3fi)^)+ I (4D 
\lb o 

+2D^ + D' + 8D a{(3) + 8{a{(3)Y - 8(/3/i)3) Nbg^^ Ubg{(3), 

L2(/3) = (^-^{16D^-12D'^-4:D^ + 3D^-{32D + 72D^-8D^-1SD^) f3fi 
+ {192D - 48^2 - 24L)3) {(3fxY - (768 - 608D - 240^2) {(3fiY - (640 

-432D) (/3/i)^ - 96 ((3^?) + — UD^ - 2D'^ - (8D + lOD^ - 13D^) f3fi 

' 16 

+ (32D - 18L>2 - 5/}=^) (/3/i)2 + 28D (/3/i)^ + 8 (/3/i)^) iVBG(/3) + ^ (8/^' 

64 

+20L'^ + 14D^ + - (16D + 56^2 + 76D^ + 26D^) - (64 - 80^ 
-1281)2 - lOOD^) (/3/i)2 + (96 - 80^ - 208^2) (/3/x)^ + (64 + 224D) (/3/i)^ 

-96(/3/i)^) (iV5G(/3))') f/BG(/3). (4.21) 

The knowledge of the explicit form of the generalized partition function fl4.18p and the 
Hill energy (14.201) enables us to compute the average number of particles and the internal 
energy. First the entropy of the system is calculated from the generalized partition func- 
tion fl4.18p through the use of (12.201) . In the next step the explicit form of the Hill energy 
fl4.20p is used in conjunction with the entropic expression in the relation = if ^ - TS 
to obtain the free energy. From the free energy the average number of particles can be 

found through the use of a specific form of (12.221) i.e., Nq'^^ = — . The explicit ex- 

pression corresponding to the average number of particles in the third constraint Ng 
is 

TVf = Nbg{(3) + (1 - g) 0Ii(/3) + (1 - qf Ol2(/3) + ■ ■ ■ , (4.22) 
where the terms 9Ti and 0^2 are 

ati(/3) = -liVBG(/3) ((2(f + 4,D)-4iVBG(/3) ((f -)) 

^2(/3) = iVBG(/3) (12 ()' + 80 0=^ + 12(8 -7D)(/3/i)2- 6/^2(2 _/))2^^ 

96 - - 

+D{1 -DK4 + D)) - i {Nscm^ (4 {f + 8 (f + 8D{f 

-16 ,D(4 + D)) + ^ {NBom' (24()' - 16 /3/x(l + /3/iX 

16 

-8(1 - 8D){(5iif + 8(2 - 3D'^)l3iJi - D{A -2D- 5D^)). (4.23) 

The calculated results corresponding to the Hill energy fl4.20p and the average number of 
particles f l4.22p is used in the relation L^^'' = f/^^^ — jj^Nq^^ to find the internal energy. The 
perturbative series corresponding to the internal energy is 

f/f (/3) = UbgW) + (1 - g) Wi + (1 - g)2 W2 + ■ ■ ■ , (4.24) 
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where the perturbative terms are hsted below 

Ui{(5) = ~{{f3fif -2f3fi + D{2 + D) + 2{D+2 f3fi 

-2{1 + D))Nbg{I3))Ubg{/3), 

U2{(5) = ^{^{I2{f + AS{l3iif -lQ{l3iif + mD{l- D){(3iif + 12D{2- D^)(5ii 

+D(4 - 3D - D^)) - 4 Nbg{P) (24 {P^f + A2D {P^if - SSD^ 

-QD{1 -D){2 + D) -_(36 [P^if + 12D [P^if + 3(16 + D^) 

+3 (iVBG(/3))' (32 (f + 4(8 + D^)J3^i + 16 {I3^if + 32(5 + 2D) {(5^if 

-(32 + 32L) - 56^2 - 2D^)(3ii + L> (8 + 20 D + lAD"^ + D^)) f/BG(/3)(4.25) 

The average number and the internal energy can also be calculated in a different manner 
wherein the basic definition of these quantities in fl2.llf2.12l) is used to find the implicit 
expressions. The implicit expressions are then solved to arrive at the final explicit forms. 
The results obtained through this method agrees precisely with those obtained with those 
in (14.221) and fl4.24p . This is a nontrivial consistency check on the results obtained for the 
physical quantities. 

The entire calculation has been carried out directly in the third constraint picture and 
consistently upto second order in the expansion parameter (1 — g). However the process 
is not limited only to the second order and it can be extended to any arbitrary order in 
the expansion parameter. 

B. Extreme relativistic classical ideal gas 

The perturbative study carried out for the nonrelativistic classical ideal gas can also be 
extended to the system of extreme relativistic classical ideal gas. Since the calculational 
procedure is same as in the nonrelativistic case detailed in the previous section, we present 
only the final expressions corresponding to the internal energy and the average number 
of particles. The average number of particles is 

(/3) = Nbg(/3) + (1 - g) Ni(/3) + (1 - qf MP) + ■ ■ ■ , (4.26) 
where the term Y^bg = exp(/3/i) V A and the coefficients of perturbation are 

MP) = -l^BG{P){{Pf^)' + 2{l + D)f3^,-D{l-D)) + {}iBGm'{{P^^)' 
-{2 + D) (3IJ + 3D-1), 

MP) = Nbg(/3) (3 iPfi)' + (20 - 12D) {f3fif + (24 - 42D)(/3/x)2 

-12D(1 -2D- D^) (511 + 2D- 3D^ - 2D^ + 3D^ - ^ i^BciP))^ ((/3/i)^ 
-2D (Pjj,)'^ - (8 - lOD) (/3/i)2 - (6 - 12D + IQD^ - 2D^) 
-D\l -6D + D')) + ^ (Nbg(/3))' (3 - (9 - 6D) (/3^)3 (4.27) 

-(2 - 28D - 3D^) {/3fi)^ + (22 - 20D - 19D^) (3fi + 2- IID + 20D^ + D^). 
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The final expression for tlie internal energy is 

Uf = Ubg(/3) + (1 - g) Ui(/3) + (1 - qf MP) + • • • , (4.28) 
and the coefficients of perturbation are listed below 

Ui(/3) = -{{(5^if-2D^3^x + D{l + D))\JBG{P) + ^BG{PW^^f-{^ + D)^3^ 
+(1 + 3D))Ubg(/3), 

MP) = (^(3(/3/i)^ + (8-12D)(/3/if -6D(l + 3D)(/3/i)2-12D2(D + l)/3/i 

+2D + + lOD^ + 3D^) - i Nbg(/3) - 2(1 + (/3/i)' 

-(4 - 2D) (/3/i)2 - (4D + 14^)2 - 2D3) (3fi + 2D + 7D^ - AD^ + 2L)^) 

+ \ (Nbg(/3))' (3 (M' + (13 - 6D) (/3^)=^ + (13 - 32/} + 3/^2) (^^)2 

-(1 + 40L> + IQD^) + 1 + 9D + 21^)2 + D^)^ Ubg(/3). (4.29) 

and the Boltzmann-Gibbs factor is U^g = ^ N^g- The Boltzmann-Gibbs values for both 
the quantities are recovered in the g — )■ 1 limit. 

The thermodynamic results corresponding to a one dimensional model of an interact- 
ing gas known as Tonks gas described by the Hamiltonian (I3.36P can be read from the 
results of the classical nonrelativistic ideal gas. For this purpose we employ a special 
kind of grandcanonical ensemble where the free length or the total length available to the 
molecules is used instead of volume. Also the number of dimensions is set to be 1 and 
the factor M.i the one dimensional analog of M. is used. 

V Generalized ensemble 

The generalized ensemble describes a completely open system which can exchange the 
internal energy, the volume and the number of particles. The heat function is the R 
function described by R = f/ + PV — fiN and the generalized free energy is the function 
S = U -TS + PV - fiN. 

A. Nonrelativistic ideal gas 

The classical nonrelativistic ideal gas is analyzed in the generalized ensemble picture. 
A perturbative process based on the idea of disentangling of g-exponential [22] is used. 
Though the thermodynamic quantities are evaluated uniformly up to (1 — g)^ order, the 
results are presented only up to the first order in (1 — g) since some of the quantities are 
too cumbrous to be displayed up to second order. The generalized partition function in 
the third constraint picture is 

if^^ P) = E NVhUW I d""^ cl-p exp, (e + - - Rf) ) ' 

(5.1) 
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where the factor j3 = and c'-^^ is the third constraint sum of g-weights in the generahzed 

ensemble. The heat function R^^ = U^^^ + PV}^^ -TS - i^N^^^ where t/f \ and ivf ^ 
are the internal energy, average volume and average number of particles respectively in the 
third constraint. The quantities P and fi denote the pressure and the chemical potential 
of the system. 

In (15. ip the phase space integrals can be carried out after the due substitution of the 
perturbative expansion of the g-exponential in terms of the ordinary exponential. Since 
the volume states are very close to each other, we approximate the volume summation 
by an integral. Employing the correct counting of volume states advocated in [25I2S], we 
carry out the volume integral. Finally we assume that the summation over the number 
of particles is convergent atleast, in the neighbourhood of g = 1 and, thus arrive at the 
final perturbative series of the generalized partition function 

Tf)(/i,P,/3) = Pi(rf^) exp(7^) Tbg0), (5.2) 

where IZ = (3 R^q^ and T>i{dp) is the derivative series defined in (14. 4 p . The partition 
function and the average number of particles in the Boltzmann-Gibbs framework is 

— +1 
(5.3) 

where the factor uj{P) is considered to be less than 1 so as to maintain the positivity 
of the partition function. From (15. 2 p the perturbative series corresponding to the third 
constraint generalized partition function upto {1 — q) order is computed 

Tf(/x,P,/5) = exp(7^)TBG(/3)(l + (l-g)^((2(a(/3))2 + 4P/3a(/3)-4(P + ^)/3 
+ (4 + ?,D))Nbg0) + (2 {a{(3)f - Aa{(3) + 2) {NBcm?) +■•■)• (5-4) 



Applying a similar perturbative technique, the sum of g-weights in the generalized ensem- 
ble picture reads: 

c(3)(;x,P,/3) = ^ -q V^idp) exp(7^) Tbg0), (5.5) 

TS')(/i,P,/3) 



where the derivative series 'D2{di3) is defined in (l4.1Up . Employing the derivative series, 
the perturbative expression of the sum of g- weights upto first order in (1 — g) is 



c(3)(/3) = ! + (!- g) (((a(/3) + l)-TBG(/3))c^(/3)-TsG(/3)) InT bg0) + ■ ■ ■ ■ (5.6) 

Adopting a procedure identical to the one used for the generalized partition function and 
the sum of g-weights, the Rq heat function can be written as: 

Rf(/i,P,/3) = ---i V,{dp)((^ + l\ i-^) exp(7^) iV5G(/3) Tsg(/3). 

(t(^)(/.,P,/3)) VV2 )p J 
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The perturbative series of the Rq upto (1 — g) order found from (15 .yp is 

Rf (/i,P,/3)= ((^ + 1^ j-fi^ NBG0) + {l-q) (Ri + R2 + R3) + ---. (5.8) 
The perturbative coefficients of (1 — g) term in fl5.8p is 

Ri = (^i?'/3(l-a(/3)) + ^((4(D + 2)(D + 4)p-i/3-2-(Z} + 2)((D + 4)(D + 6) 

+8 (i? + 2 fi)p-^)^-^ + {R + fi) fiP-^ + 3 /i - 4 ((D + 2) (i? + 4 
-{2R^ -6fx-3D + 6D{D + 6) fi + 8{R + fi) 

+UJ0) (2 {D^ + 40D^ + 32D + 32 + 2{D + 2)R^ + A{D + 2){R + 2 fi)) P~^p-^ 
-4 {D + 2){D + 4) _ 4(2(P)2 + 5/) + 3) ^ + p + 2)(3P' + 10) 

+4 {R + /i) ^P"^) + 8((2P' + 4) P + 3P> /i) - ((P + 2iif 

R2 = ^(^{D + 2f r^ + {D + 2){{D + 2){3D + 8) + 12 fi^)- 2 {12 fi^P^ + 4 fi^p-^ 

+4(9P'2 + 34P) + 32) /iP2 + 2(P) + 2f RP^) (3^ + 4 {{9D + M)/!^ 
+2{D + 2) RfiP^)l3^ - 12(P) + 2) V) - ii^D^ + 20P'2 + 44P) + 32) P^^ 
-2 ((9P>2 + SAD + 32) + 2(P) + 2f P) /3P^ - 4 (4(P) + 2)P 

+(9D + 14) /i) P2/3V + (2P + 3/i) /3Vp2))Tbg(/3) (iVscl/S))', 

R3 = l(3uj0){{D + 2fP'^(^^-6{D + 2f^iP^f3^ + 12{D + 2)fi^f3^P^-8fi^(3^P'^) 
8 \ 

-3{{D + 2)3 /32p3 - 6(P' + 2) + 12(P' + 2) /3V^^ 

-8 /3 V^')) Tsg(/3) iNBG0) f. (5.9) 

The generahzed partition function fl5.4p . the sum of g- weights fl5.6p . and the Rg^^ heat 
function (15. 8p . form a set of coupled imphcit perturbative expansions. The fact that these 
equations are uncoupled at the lowest perturbative order allows us to obtain an explicit 
form of these quantities through the use of a recursive procedure. The final explicit form 
of the generalized partition function computed upto first order in the expansion parameter 
(1 - q) is 

Tf )(/3) = Tbg(/3) exp(/3P) (1 + (1 - g)(Ti(/3) + T,iP) + T,iP) + T,m + ■■■), 

(5.10) 
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and the coefficients of perturbation are listed below: 

T^{/3) = -^[A{2{a{/3)r-Aa{/3) + D)-{{2 + D)^{{2 + D){A + D-/3f,) 

-4(/3/i)2) - 4(4 + QD + D^ + 2(2 + D) 13^) In Tbg(/3)) Tbg(/3) 
+2 {D{2 + D) - 4 {(5^f - 4 (2(2 + D) - 2(1 - D) (5^ 

-(/3//)2) \uT BoifiW BG{fi)f + 16 ^^(1 + aim^BGmf) NbgW), 

T2(/3) = -^((6 + 3D-(8 + D)/3/i-(/3^)2-2(l + i?-/3^)lnTBG(/3)) 

+ (16 + 8D- 2D2 - 16/3P - (8 - lOD + + 8/3P - 4D pP)Pn 
-(8 + 4P) + 8/3P) (/3/i)2 + 12 (/3/i)3 + 4(6 + - 2(/3P)-^ - (4 - P>) 
+ (/3/x)2) In Tscm Tbg(/3) - 4 ((/^P)-^ - 4 (/3P)-^a(/3) 

-(a(/3))2) {Tsam') {NsGmr, 

T3(/3) = -2^^^1nTBG(/3)-^(8 + 6D + D^-(16 + 4D + 4/3P)/3^ 

+12 (/3//)2 + 4(1- a(/3)) In TbgW) + 4(«(/9))' + (/3P)-'a(/3)) (A^sG(/3))^ 
T4(^) = -^(a(/3))^(7V^G(^))^ (5.11) 
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Proceeding similarly the final explicit form of the heat function Rg is obtained upto 
(1 — g) order 

+ ^"'^) + ^^~^^ (Ri+R2 + R3) + --- , (5.12) 
where the perturbative coefficients may be listed as 
Ri = ^((P» + 2)(P> + 4)r'-p + 2)//-/3//')(TBG(^))'lnTBG(^)Ar5G(^), 

Y^(^{iD + 2){3D + 10)P -{D + 2)(3/x - (D^ - IQD - 40) P) P/3 

+2 (24 At - (P> + 2)(3P) - 16) P) /iP^^ + 4(3P) - 20) /J^P^p^ - 8 //^P^^^ 
+2 (8 - 3(P> + 2)(P> + 4) P) ^/x + 12(P' + 2) P^iif - 8 P(^/i)=' 
-{D + 2)(16/x + (D + 4)(P» + 6) P)) Tsg(/3) - 2((P' + 2)^(P> + 4)^' 
-2(P) + 2)(3P' + 8) /i^ + 12 (P) + 2) I3i2^ - 8 /3V + 2 ((D + 2)2/3"! 

-4(D + 2)// + 4/3//2)lnTBG(/3)) (Tbg(/3))') (iVBG(/3))^ 

^ (2 (P) + 2)((9P) + 10)/x + 2(P)2 + 3D + l)P)P^f3^ - D{D + 2)(3P) + 8) P^/? 
8 \ 

-4 ((gP) + 14) + (Pi + 2)(3P + 4) p2) ;^/33 + 8 (3// + 2(3P + 5) P) /x^P^/^^ 
-2PV/3')(A^bg(^))'. (5.13) 
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To compute the other thermodynamic quantities hke the internal energy, average volume 
and the average number of particles we need to calculate the free energy via the expression 
Sq'^^ = ujf^ + PVq^^ — TS — fiNjf^ = Rg^'' — TS . The entropy which can be read from 
the generalized partition function f lS.lOp using the relation fl2.20p . and the perturbative 
expansion of the heat function Rq are used in this process. From the free energy, the 

(3) dS^^^ 

average volume can be found through the use of the thermodynamic relation Vg = -^p-- 
The perturbative expansion of the average volume upto (1 — g) order is 

Vg^'^ = + ^ (2(a(/3))'-2«(/3) + 4 + D + 3(2(a(/3))2-4a(/3) 

+8 + D) Tbg{(3) + 4 {{com' + {a{(3) - 3) c^(/3) - - 2)) In Tbg(/3) 

(1 + um Ham' - + 2))yBG(/3) + ■ ■ ■ , (5.14) 

where Vbg{P) is the average volume in the Boltzmann-Gibbs formalism. Similarly the 

(3) d£^^^ 

average number of particles can be arrived through the use of the relation Ng = 
The final expression of the average number as a perturbative series upto first order in 
(1 - q) is 

iVf = iV5G(/3) + (l-g) ^ ((4(/3/i)' + 8/3P«(/3)-D(2 + D) + 8(«(/3) (5.15) 
+1) In Tbg(/3)) iVBG(/3) - (12 (/3/i)' - 4D - 8 /3P(a(/3) + 1) - D(2 + /}) 
+8 (a(/3) - 1) In Tbg(/3)) {^30^? + 4 (2 - /3P - 1) [^30^?] + ■ ■ ■ • 



Using the perturbative expansions of the heat function fl5.12p . the average volume (15.141) 
and the average number of particles fl5.15p . the internal energy can be computed via the 
legendre transformation Rg'^'' = uf"^ + PVg^^ — fiN!f\ The computed perturbative series 
of the internal energy upto (1 — g) order reads: 

[/f = UBGi(3)-{l~q)lrBG{(3){i8lnTBG{(3) + a{(3)-l)u{(3) + 8{a{(3) 

-2) In TbgW) + 4 (/3^ - 2)a(/3) + 2Df3fi-{D + if) Ubg{P) + ■■■ ■ (5.16) 

The quantities evaluated above can also be arrived at using the basic definition of the 
internal energy, the average volume and the average number of particles. For this pur- 
pose we employ the basic definition of these quantities to find the corresponding implicit 
equations which are then solved recursively to obtain the final explicit expressions. The 
results obtained through this method perfectly agrees with the results displayed above. 
All the quantities recover the corresponding Boltzmann-Gibbs results in the g — )■ 1 limit. 

B. Extreme relativistic ideal gas 

The extreme relativistic ideal gas has been studied using the perturbative formalism. 
The actual process of calculation is same as in the nonrelativistic case and so we present 
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the final results directly. The perturbative expansion corresponding to the the average 
number of particles is 

Nf ) = ^BG - Tsci^) [d + D'- 2D /3/i + 

-2{l + D- f3fi) InTBGiP)) + ■ • • • (5.17) 

The factors TbgW) Nbg{(3), are the partition function and average number in the 
Boltzmann-Gibbs statistics respectively and their expressions are given below 

1 £(/3) 1 

W« = Y3«]fy. ""--r^y «« = -p(MA-(-j (5.18) 

is the average number in the Boltzmann-Gibbs statistics and ^(/3) is assumed to be less 
than 1. The average volume in nonextensive statistical mechanics is 

= VBG-^-^-^rBGm{{2 + 3D + D'-2{l + D)Pf,-{(3fif (5.19) 
-2(2 + D-Pfi) lnTBG{P))'^BG{P) -2{1 + D~Pfi- \nTBG){^BGm') +■■■ ■ 

where the Boltzmann-Gibbs value corresponding to it is Vbg = jp ^BGi(3). A similar 
perturbative series corresponding to the internal energy is 

= \JBG-^-^^rBGm{{2 + 3D + D'-{l + D)f3f,-{^fiY (5.20) 
-2(2 + D-Pfi) lnTBG{l3))VBG{P) - 2(1 + D - - lnTBG){^BGm') +■■■ ■ 

The Boltzmann-Gibbs internal energy used in (15:^ is VbgW) = j ^bgW)- From the 
results displayed above we can observe that in the g — 1 limit we recover the standard 
Boltzmann-Gibbs values. 

The Tonks gas described by the Hamiltonian fl3.36p can also be studied through the 
generalized ensemble. The results corresponding to such a model of interacting gas in the 
free length generalized ensemble can be obtained by substituting D = 1 in the results 
obtained for the nonrelativistic classical ideal gas. 

VI Adiabatic ensemble 

The adiabatic class of ensembles is used to describe a system of particles for which the 
thermal equilibration is with respect to the heat function. In the adiabatic class the 
individual members of the ensemble have the same value of the heat function though they 
can be at different temperatures. There are four different adiabatic ensembles namely 
microcanonical (A^, V, U), isoenthalpic-isobaric ensemble {N, P, H), the ensemble with 
number fluctuations (/i, V, L), and the ensemble with both the number and volume 
fluctuations {fi, P, R). In the current section we present a unified framework to describe 
the adiabatic class of ensembles. 
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At any given instant in time the microstate of a system consisting of N particles is 
represented by a point in the 6A^ dimensional phase space, which is comprised of 3A^ posi- 
tion coordinates and the 3A^ momentum coordinates. Since the position and momentum 
coordinates evolve with time, the representative point moves in the phase space. The 
motion of the representative point traces a trajectory of constant heat function in the 
phase space. The points which lie on the surface of the constant heat function curve are 
the various microstates corresponding to the macrostate of constant heat function. As 
the number of microstates is very high and lie very close to each other, the surface area 
of the constant heat function curve can be considered as a measure of the total number 
of microstates. The surface area of the constant heat function curve in an adiabatic 
ensemble is 

Iff ^ 
(](Xi,X2,io) = J] ^— ^ / / 6{n + J2^iXe-^)U^''r,d^P^■ (6-1) 

The computation of area of the constant heat function is difficult, and so we usually 
calculate the volume enclosed by such curve. The phase space volume enclosed by the 
constant heat function Sj curve is 

N 

X2, ^) = Y1 ]^nli5^ / / + E ^^^^ - ^) n d""^^ "^"P^- (6.2) 

Xiey ■ -^P' {£} i=l 

The surface area of the heat function curve (16. ip and the volume (16. 2p enclosed by it are 
related via the relation 

r(Xi,X2,i?) = ^fi(Xi,X2,i3). (6.3) 

Based on the kind of adiabatic confinement there are four different definitions of the 
entropy. A unified definition of the entropy for an adiabatic nonextensive system is 

^,(Xi,X2,io) = k ln,r(Xi,X2,i^), (6.4) 

where the i7(Xi, X2,io) is the volume enclosed by the curve of constant heat function and 
a measure of the number of microstates. The temperature of a general adiabatic ensemble 
is defined via the relation 

(dSA-' {E{X^,X„m' .... 

In the usual procedure of calculation the heat function can be computed through the use 
of (16.51) and the relations (16. ip and (16. 2p . Using the heat function the specific heat can 
calculated through the expression 



a 



'"^w dT 



(6.6) 



The extensive thermodynamic variable whose intensive counterparts are held fixed can be 
obtained from the relation 

1 d 

X£;g(Xi,X2,io) = -- —Sg{Xi,X2,S)). (6.7) 
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The volume and the average number of particles can be calculated in the {N, P, H) 
and (/i, V, L) respectively through the use of (16.71) . In the case of (/i, P, R) ensemble, 
both the volume and the average number of particles can be calculated using (16. 7p The 
expectation value of an observable O for any adiabatic ensemble can be defined in a unified 
way through the relation 

TV 

- mk^) LL" '^''^V, - ^) H ^'r, d-p.. (6.8) 

The internal energy is the natural heat function of the microcanonical ensemble and so 
it can be directly evaluated, whereas in the other three ensembles, equation (16. 8p is used 
to calculate the internal energy. From the internal energy the specific heat at constant 
volume can be evaluated using 

dU, 



C, 



1\V 



dT 



(6.9) 

V 



A nontrivial way of checking the internal energy can be accomplished by using the Leg- 
endre transform (11.31) to find the internal energy. 

A q generalization of the equipartition theorem is developed in the unified framework. 
Let t/i be a phase space variable which can denote either the position or the momentum 
Pi {i = 1, 3A^) coordinate. Calculating the expectation value of yi |^ we arrive at 

Vi —J = kT (r(Xi,X2,io))^-' 5,,. (6.10) 

Equation (I6.10p obtained above is the generalized form of the g-nonextensive equipartition 
theorem. It is interesting to note that the nonextensive generalization of the equiparti- 
tion theorem depends on the factor Z'(Xi,X2,i3) which is a measure of the number of 
microstates corresponding to a given macrostate. However we notice that in the extensive 
g — > 1 limit the dependence on the number of microstates vanishes. The form of (I6.10p 
suggests that the the expectation value described has a nonzero value only when i = j. 
When the phase space variable i/i is set to be the coordinate rj then we get a specific form 
of the equipartition theorem 

n^) = -{n p^) = {n F,) = kT (r(Xi,X2,-fi))^-' , (6.11) 



where we have used Hamilton's equation of motion |^ = —pi. Following an identical 
approach for the situation where the generalized variable is set to be the momentum 
variable Ui = pi we arrive at 

81-1 \ 

Pi ) = {Pi q^) = kT (r(Xi,X2,-fl))''"' . (6.12) 

We observe that (I6.12p is nothing but twice the expectation value of kinetic energy. From 
(16. lip and (I6.12p it can be proved that systems whose Hamiltonians can be cast in the 
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form "H = Pi + J2i Ql through a canonical transformation the expectation value 

of the Hamiltonian is 

n = ^kT (r(Xi,X2,i3)r\ (6.13) 
The virial theorem in the nonextensive framework obtained from (16.111) reads: 

/ \ 

{Y^q^p^) = -^N kT {S{X,,X2.S))y-\ (6.14) 

^ i ' 

The unified framework derived above is made use of in the following sections to demon- 
strate the isoenthalpic-isobaric ensemble {N,P,H), the ensemble with number fluctuations 
(/i,V^,L) and the ensemble with both number and volume fluctuations (/i,P,R). Though 
all the relations in the microcanonical ensemble {N,V,U) can be read from the equations 
in the unified framework above, we do not discuss the {N,V,U) ensemble, since a detailed 
study has already been carried out in [2T] . 



VII Isoenthalpic-isobaric ensemble 

The study of adiabatically confined systems with constant enthalpy is done using the 
Isoenthalpic-isobaric ensemble. In the following discussion we study the classical nonrel- 
ativistic and extreme relativistic ideal gas. 



A. Nonrelativistic ideal gas in (A^^, P, H) ensemble 

The nonrelativistic ideal gas is studied in the isoenthalpic-isobaric ensemble. The volume 
enclosed by the surface of constant enthalpy is calculated in the following manner. First 
the Hamiltonian corresponding to the nonrelativisitic gas (13. ip is substituted in equation 
(16. 2p and the resultant expression reads: 

^(^. ^. ^) = E ]^ /_ { e (^E I; + - ^) n ""p.- (") 

The momentum integration in (17. ip is geometrically the volume of DN dimensional sphere 
of radius \/^2/m{W^^lW) and is equal to 

„.I, = d°>. = (H - PV)^. (7.2) 

Ji:?''^<H-pv i (— + 1) 

Substituting the result of the momentum integration (17.21) in (17. ip the expression for the 
volume is 

In the next step the integral over the position coordinates is carried out in (17. 3p and this 
yields 

g(iV.P.g) = :^ ^^ Y^V^(H-PV)'^. (7.4) 



V 
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Finally the summation over the volume eigenstates is considered. Since the volume states 
are very closely spaced the summation is approximated by an integral. But as we have 
discussed earlier an integration leads to overcounting of volume states. To overcome this 
we employ the shell particle method of counting of volume states proposed in [25|26] . 
wherein only the distinct minimum volume states needed to confine the particle is taken 
into account. The volume enclosed by the isoenthalpic curve in phase space is 



E{N,P,H) = M'' (- 



1 



N 



r(Di 



1^ 



The surface area of the constant enthalpy curve in the phase space is 



n{N,P,H) = M'' {- 



r(Di 



(7.5) 



(7.6) 



From the volume of the isoenthalpic curve in the phase space the nonextensive entropy 
of the classical ideal gas in this ensemble is 



S,{N,P,H) = k 



{m- (^) 



r(Di+i) 



v 



(7.7) 



Using the definition of temperature (16. 5p . the heat function corresponding to the ensemble, 
the enthalpy is found 



S.1 



P 



{l-q)NZi 



{7.t 



The specific heat at constant pressure evaluated from (17. 8p through the use of (16.61) is 



r(Di 



(7.9) 



The average volume calculated using (16. 7p can be rewritten to obtain the equation of state 
given below 



PV„ 



r(S)i 



(7.10) 



The internal energy evaluated through the definition of the expectation value of the ob- 
servable ( 16. 8 p is 



DN 



M 



N 



2/3 \r(2)i + l) \P 



N 



(7.11) 



The calculated value of the internal energy has also been checked through a different 
procedure using the Legendre transformation Hg = Uq + PVg. From (17. lip the specific 
heat at constant volume found using (16. 9p is 



c, 



DN 



i\v 



k £i 



M 



N 



r(2)i + i) \p 



(7.12) 



30 



Using the results corresponding to the specific heat at constant pressure (17.91) and the 
specific heat at constant volume f l7.12[) the ^-generalization of the Mayer's relation is 
obtained 



C,|,-C,|, = 7V*£. 1^^^ ^-j ^-j j , (7.13) 

The ratio between the specific heat at constant pressure (17.91) and the specific heat at 
constant volume fl7.12p 

'' = 7?r = i + ^ 

is interestingly independent of q and A^. We have observed the same result in the 
isothermal-isobaric ensemble f l3.2ip . In the extensive g — )■ 1 limit the standard Boltzmann- 
Gibbs results are recovered. 



B. Extreme relativistic ideal gas in [N^P^H) ensemble 

The extreme relativistic classical ideal gas described by the Hamiltonian fl3.23p is studied 
in the {N, P, H) ensemble. The volume enclosed by the surface of constant enthalpy is 
given by 



N 

V ' "^'"^ -^P' i i=l 



(7.15) 



In (I7.15P the result of the momentum integration is 



T 

er-'-p 



Y:?''v^<H-pv 



r {DN + 1) 



27rT T{D) 



N 



{H - PV) 



DN 



(7.16) 



The expression (I7.16P is substituted in (I7.15P and the resulting expression is integrated 
over the position coordinates. Replacing the summation over the volume eigenstates by 
an integration to consider the continuum nature of volume and adopting the shell particle 
method of counting of volume states the final expression for E{N, P, H) is 



r(iV, P,H) = A 



N 



N 



H^2 



(7.17) 



pj r(D2 + i)' 

which enables us to calculate the surface of the isoenthalpic curve via the relation (16.3 



n(N, P,H)=A 



N 



r(S)s 



(7.18) 



The nonextensive entropy (16.41) of the extreme relativistic ideal gas obtained from (I7.17P 
is 



S,iN,P,H) = k 



(a- (1) 



N 11^2 



r(2)2+i) 



1-9 



l\ 



l-q 



(7.19) 
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Using the definition of temperature in tlie isoentlialpic-isobaric ensemble f l6.5p in conjunc- 
tion witli tlie expression for tlie entropy (17.191) we find tfie entlialpy 



V(r(D2))^-V \P) 



(7.20) 



Tlie specific lieat at constant pressure obtained from tlie enthalpy via (16. 6p is 

(i) ) . (7.21) 

The computed expression for the average volume using (16. 7p is rewritten to give the 
equation of state 




{T)2f' fiY /l^^^ 



'2 



Making use of the formula for the expectation value (16. Sp . the calculated form of the 
internal energy is 

DN / A" n x" 
- - (r(57TT) (?) is) ) ' f^-^^' 

and this perfectly agrees with the result found using the Legendre transformation H = 
U + PV . From the internal energy the specific heat at constant volume can be found: 

The ratio between the heat capacities at constant pressure and constant volume 

is found to be independent of both the nonextensive parameter and the number of particles 
A^, similar to the result observed in (I3.35p . The enthalpy, internal energy and the specific 
heats go to their respective Boltzmann-Gibbs value in the g — >■ 1 limit. 

The thermodynamic quantities corresponding to the Tonks gas (I3.36P can also be com- 
puted through a procedure identical to the one employed in the case of the nonrelativistic 
and extreme relativistic classical ideal gas. The results can be easily read from the ex- 
pressions corresponding to the nonrelativistic gas through the use of the free length Lf.q 
in the place of the volume V and setting the dimension D = 1. The equation of state of 
the Tonks gas thus obtained is 



It can be seen that equation (17.260 recovers its standard Boltzmann-Gibbs result PLj.q = 
N kT in the g — t- 1 limit. 



32 



VIII Adiabatic ensemble with number fluctuations 



The adiabatic ensemble with number fluctuations was first discussed in detail in [lOj in 
the context of Boltzmann-Gibbs ensemble. This ensemble is the adiabatic counterpart 
of the grandcanonical ensemble and the Hill energy is the corresponding heat function. 
Though the integration with respect to the position and the momentum coordinates can 
be carried easily, it is not possible to carry the sum over all possible values of A^. So in 
this section we give only a formal expression corresponding to the phasespace volume, 
surface area, the Hill energy and the equation of state. 



A. Nonrelativistic Ideal gas 

The volume of phasespace enclosed by the curve of constant Hill energy L after carrying 
out the integrations over the position and the momentum coordinates is 

^(^' ^' L) = E M^v(I2N^,^ (L + Z^^)"^- (8-1) 

The surface area corresponding to the phasespace curve of constant Hill energy can arrived 
at through the use of (16. 3p . The equation of state describing the classical ideal gas in this 
ensemble is 

= f E^ jv,rL + il (L + A'V)™ (8.2) 

N=0 ' V 2 / 

The formal expression for the average number of particles arrived at through the use of 
(EZD is 

§ IlfmYjm) ^" (L + fiv)'*^-'. (8.3) 

B. Extreme relativistic Ideal gas 

For an extreme relativistic classical ideal gas described by the hamiltonian fl3.23p . the 
phase space volume calculated is 

The equation of state and the average number of particles of the extreme relativistic gas 
corresponding to this adiabatic ensemble are 

kT ^ 

PV^ = — Yn——^ (L + /iiV)^^, (8.5) 

OO 



- g n^uL-.) (L + .iV)-'. (8,6) 
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In the g — )■ 1 the above quantities reduce to the standard Bohzmann-Gibbs equations. 
Substituting D = 1 and replacing the volume by the free length in the thermodynamic 
expression for the nonrelativistic ideal gas, we can arrive at relations corresponding to the 
Tonks gas. 

IX Adiabatic ensemble with volume and 
number fluctuations 

An adiabatically confined system which can exchange both volume and particles with 
the bath is described by this ensemble. In the classical Boltzmann-Gibbs case such an 
ensemble was first described and discussed in detail in ^llj. The heat function of the 
system is the R heat function and this is the adiabatic counterpart of the generalized 
ensemble. Similar to the previous ensemble, the N summation is difficult to carry and so 
we present the results formal sum. 



A. Nonrelativistic Ideal gas 

In the ensemble with both volume and number fluctuations, the phasespace volume of the 
classical nonrelativistic ideal gas enclosed by a curve of constant R is 

°° / 1 \ ^ 

r(/.,p,R) = X; [p) (R + ^ivf^ (9-1) 



Using the expression for the phase space volume fl9.1|l . the equation of state and the 
average number of particles is obtained 



N=0 

oo . ,^ 



^. = f E^rWT (R + z^^)"^" (9-3) 



B. Extreme relativistic Ideal gas 

In the case of an extreme relativistic ideal gas, the phase space volume enclosed by a 
curve of constant R is 

which leads to the equation of state and the average number of particles given below 

kT ^ / 1 \ ^ 
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The g — )■ 1 limit is respected in all the above mentioned cases. The thermodynamic 
expressions corresponding to the Tonks gas can be arrived directly from the relations 
given for the nonrelativistic gas by substituting D = 1 and replacing volume by the free 
length. 

The summation over the number of particles could not be evaluated in both the ensem- 
ble with the number fluctuations and the ensemble with number and volume fluctuations. 
But the thermodynamic quantities in both these ensembles may be obtained through the 
use of a molecular dynamics simulation. 

X Remarks 

A comprehensive investigation on the different kinds of ensembles has been carried out 
in the present work. We notice that there are two different classes of ensembles namely 
the isothermal class and the adiabatic class. The thermal equilibration in the isother- 
mal class happens with respect to the temperature. Isothermal class comprises of the 
canonical ensemble, the isothermal-isobaric ensemble, the grandcanonical ensemble and 
the generalized ensemble. We notice that the thermodynamic relations corresponding to 
these four ensembles can be presented through a unified formulation. Such a formulation 
is evolved for both the second and the third constraint formalisms though the latter is 
the currently accepted formulation of nonextensive g-statistical mechanics. The develop- 
ment of unified formulation corresponding to the second constraint has been done with a 
view to apply nonextensive q statistics in computational simulations like Montecarlo and 
molecular dynamics. It is a well known fact that the probabilities in the third constraint 
are implicit quantities. In calculational procedures involving simulations it is easier to 
find the thermodynamic quantities in the second constraint and later transform them to 
the third constraint. Towards this end a generalization of the temperature dependent in- 
terrelation between the second and the third constraint introduced in Ref. [12] has been 
accomplished. 

The unified formulation developed is applied to study the isothermal-isobaric, grand- 
canonical, and the generalized ensembles through specific examples viz the classical non- 
relativistic, and, extreme relativistic ideal gas models. Since the canonical ensemble has 
already been studied in detail in Ref. P^fTB] it has not been considered in our present 
work. In the isothermal-isobaric ensemble, the implicit equations corresponding to the 
generalized partition function and the enthalpy are solved to obtain their corresponding 
explicit expressions. From these quantities the average volume, the internal energy, the 
specific heat at constant pressure, and, the specific heat at constant volume are obtained. 
All the above mentioned thermodynamic quantities are obtained as exact results. The 
specific heats were dependent on the nonextensivity parameter q and the temperature 
in both the models. Interestingly we notice that the ratio between the specific heat at 
constant pressure and the specific heat at constant volume is independent of both q and 
A^. A similar exact computation of the gas models in the grandcanonical and the gen- 
eralized ensembles were found to be difficult. So, we employed a perturbative procedure 
developed in [T8|23] based on the disentangling of g-exponential [22]. The perturbative 
series corresponding to the generalized partition function, the sum of g-weights and the 
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heat function are obtained up to a predetermined order in the third constraint formahsm. 
The determined imphcit quantities are then solved recursively to obtain the final explicit 
forms which are used to find the thermodynamic results. The fact that these quantities 
arc uncoupled and exactly solvable in the q = 1 limit facilitates the recursive procedure. 
In the grandcanonical ensemble, the final results are computed up to second order in the 
expansion variable {1 — q) whereas in the generalized ensemble the results are displayed 
only up to the first order. Though the results have been calculated and displayed only up 
to a particular order in the expansion parameter, the procedure can be extended to any 
arbitrary order in (1 — q). 

A similar unification procedure is used to create a generalized formulation for the 
adiabatic ensembles. The microcanonical ensemble, the isoenthalpic-isobaric ensemble, 
the adiabatic ensemble with number fiuctuations, and the adiabatic ensemble with num- 
ber and volume fiuctuations belong to this class. The nonextensive q generalization of 
the equipartition theorem and the virial theorem have been obtained in a unified sense. 
The {N,P,H) ensemble, {^,V,L) ensemble and the {fi,P,R) ensemble have been studied 
through the examples of classical nonrelativistic ideal gas and the extreme relativistic 
ideal gas. The above mentioned ideal gas models have been solved in the isoenthalpic- 
isobaric ensemble. The exact expression corresponding to the entropy, the enthalpy, the 
internal energy and the heat capacities at both constant pressure and at constant volume 
were found. In the case of {fi,V,L) and the (/i,P,R) ensemble an exact solution could 
not be arrived and so the thermodynamic quantities have been expressed as formal sums. 
An interesting observation is made from the study of the ideal gas models in the vari- 
ous ensembles. The gas models were exactly solvable in the isothermal-isobaric and the 
isoenthalpic-isobaric ensembles. In both these ensembles we notice that the ratio of the 
specific heat at constant pressure and the specific heat at constant volume was found to 
be independent of the nonextensivity parameter q, the number of particles N and the 
temperature. This is despite the fact that the individual heat capacities were dependent 
on these three factors. Meanwhile the difference between the heat capacities was found 
to be dependent on these three factors. 

The choice of an ensemble to study a problem in statistical mechanics depends on 
the physical situation. The microcanonical ensemble is used to study completely isolated 
systems, whereas the canonical ensemble is used to study closed systems which exchange 
energy with the surroundings. The grandcanonical ensemble is used to study systems 
capable of exchanging particles with its surroundings. A system like colloidal particles 
or macromolecules immersed in a solvent experiences a constant pressure and so they 
can be suitably studied using the isothermal-isobaric ensembles. Depending on whether 
a fiuid under constant pressure is isothermally or adiabatically confined, we can use the 
isothermal-isobaric or the isoenthalpic-isobaric ensemble. The generalized ensemble can 
be used to analyze systems whose size varies while the pressure is maintained constant. 
Typical examples include molecular clusters in a one component gas and a system of 
bound ions or molecules on a protein molecule. Also there are situations in nature where 
an external electric field, magnetic field or gravitational field acts on a thermodynamic 
system. The construction of a class of ensembles comprising of an external field along 
with the thermal, mechanical and chemical variables, and their consequent application to 
study problems like rotating black holes and charged polymers like DNA molecules should 
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be worth pursuing. 

The unification procedure developed in the current work for the nonextensive g-entropy 
(11 .ip can be extended to other kinds of nonextensive entropies like k entropy [28] and the 
two parameter entropies [29]l30] . Such an effort will help us to understand the change in 
the structure of thermostatistics due to nonextensivity. It has been noticed in Ref. [TU] 
that there is an inherent pairing between the isothermal and the adiabatic ensembles that 
is each isothermal ensemble has an adiabatic counterpart. These pairs can be identified 
through the fact that they have the same heat function. An interrelation between the 
isothermal and adiabatic ensemble can be established through a Laplace transformation. 
The canonical ensemble and the microcanonical ensemble have been interrelated in the 
nonextensive g-statistics through a g-generalization of the Laplace transform [31j. Cur- 
rently we are working on developing a unified formulation of the g-Laplace transform to 
connect any of the isothermal ensemble to its adiabatic counterpart. The results will 
appear elsewhere. 
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